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Abstract: The study is concerned with a different perspective which the numerical solution of the singularly 
perturbed nonlinear boundary value problem with integral boundary condition using finite difference method on 
Bakhvalov mesh. So, we show some properties of the exact solution. We establish uniformly convergent finite 
difference scheme on Bakhvalov mesh. The error analysis for the difference scheme is performed. The numerical 
experiment implies that the method is the first order convergent in the discrete maximum norm, independently of 
£- singular perturbation parameter with effective and efficient iterative algorithm. The numerical results are 
shown in table and graphs. 
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1. Introduction 


This paper focused on the first-order nonlinear singularly perturbed differential equation. 

eu+f{x,u) = 0, x e (0,L], L> 0 (1) 


subject to the integral boundary condition 

u( 0) = /?u(L) + f b(s)u(s)ds + d 
“'0 


( 2 ) 


where £ is a small perturbation parameter (0 < £ « 1). a, fi, d are constant, which are independent of £. We 
assume that b(x), f(x, u) are sufficiently smooth functions in the intervals (0, L] and (0 ,L\xR. Also, 
j- > a > 0. The solution of (1) has in general a boundary layer with 0(£) near to x = 0. The first time, 

nonlocal problems were studied by Bitsadze and Samarskii [16]. Differential equations with conditions which 
connect the values of the unknown solution at the boundary with values in the interior are known as nonlocal 
boundary value problems [16]. Existence and uniqueness of nonlocal problems can be seen in [1-3]. Singular 
perturbation problems can also arise in fluid mechanics, quantum mechanics, hydromechanical problems, 
chemical-reactor theory, control theory, oceanography, meteorology, electrical networks and other physical 
models. The subject of nonlocal boundary value problems for singularly perturbed differential equations has 
been studied by many authors. For example, Cziegis [4] studied the numerical solution of singularly perturbed 
nonlocal problem. Cziegis [5] analyzed the difference schemes for problems with nonlocal conditions. 
Amiraliyev and (jakir [8] applied the difference method on a Shishkin mesh to the Singularly perturbed three- 
point boundary value problem. Amiraliyev and (jakir [9] studied numerical solution of the singularly perturbed 
problem with nonlocal boundary condition, (jakir and Arslan [14] analyzed a numerical method for nonlinear 
singularly perturbed multi-point boundary value problem. Adzic and Ovcin [11] studied nonlinear spp with 
nonlocal boundary conditions and spectral approximation. Amiraliyev, Amiraliyeva and Kudu [12] applied a 
numerical treatment for singularly perturbed differential equations with integral boundary condition. Herceg [6] 
studied the numerical solution of a singularly perturbed nonlocal problem. Herceg [7] researched solving a 
nonlocal singularly perturbed problem by splines in tension, (jakir [10] studied uniform second-order difference 
method for a singularly perturbed three-point boundary value problem. Geng [15] applied a numerical algorithm 
for nonlinear multi-point boundary value problems, (jakir [13] obtained a numerical study on the difference 
solution of singularly perturbed semilinear problem with integral boundary condition. It is well known that 
standard discretization methods for solving nonlocal singular perturbation problem are unstable and these don’t 
give accurate results for £. Therefore, it is very important to find suitable numerical methods to these problem. 
So, we use finite difference method in this paper. According to this method, we will analyze that this method for 
the numerical solution of the nonlocal problem (l)-(2) is uniformly convergent of first order on Bakhvalov 
mesh, in discrete maximum norm, indepedently of singular perturbation parameter £. In section 2, some 
properties of the exact solution of problem described in fl)-(2) is investigated. According to the perturbation 
parameter, by the method of integral identities with the use exponential basisfunctions and interpolating 
quadrature rules with the weight and remainder terms in integral form uniformly convergent finite difference 
scheme on Bakhvalov mesh is established in section 3. The error analysis for the difference scheme is 
performed in section 4. In section5, numerical examples are presented to find the solution of approximation. In 
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this paper, we use Il 5 '(x)|| 00 for the continuous maximum norm on the corresponding interval. g(x) is any 
continuous function. 


2. Properties of the Exact Solution 

The exact solution of (l)-(2) has some properties, which are needed in later sections for analysis of the 
numerical solution. Thus, we describe these properties with Lemma 2.1. 

Lemma 2.1. If > a > 0 and/(x,u) e (0 ,L]xR, the estimates 

II u(x) || < C 0 (3) 

where 

—at / —aT \ 

C 0 = |u(0)|e~ + (l - e^Ja^llFlL, F(x) = -/(x,0) 

and 

f 1 ~ ax 1 

Hit (x) || <c[l+-e—j (4) 

hold for the solution u(x), where C 0 and C are constant independent of e and h. 


Proof. We use intermediate value theorem for /(x, u) in Equation (1) 
f{x, u ) — fix, 0) dfix, u) 


and 


du 


a(x) 


u = yu, 0 < y < 1 


dfix, u) 
du 


Thus, we obtain the following linear equation, 

su + a(x)u(x) = Fix ) 

We can write the solution of the Equation (5) as follows: 


u(x) = u(0)e^ J o aCr)dT + - [ Fif)e J f ( a(r)dT 

£ Jn 


df 


In Equation (5), x — T and x = s, we obtain Eq (7)-Eq(9) 


df 


and 


u(0) = 


u(T) = u( 0 )eT^«* + - I F^ e ’(“W* 

£ J o 

u(s) = u(0)eT‘ J +- [ Ftf)e^ S ( a(0d< df 

£ J o 

d + F(f)e^ J ' a(T)dT df + -Jl bis ) [f FiOe^ aMdr df]ds 


1 - “W* - J o r bis)e ^ S ° a(T)dT ds 

Now, we evaluate to Equation (9) for a(r) > a > 0 and b* = max\h (x) |. 
First, we prove the estimate the following inequality. 


l _ p e ~ril a(T)dT 


f 

Jo 


-aT 


—1 rs ( , -aT I ‘ —as —at / 

-j bis)e~ J o aMdr ds>l-/Se~-b* j e~ds > 1 - /Se~ - b* [ 


■f 

Jo 


(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


1 — e 




> 0 

Second, we prove the estimate the following inequality. 
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d + - f F(Oe^ f t a(T)dT df + - [ b(s ) [ [ F(Oe^ a(r)dr d(]ds 
£ Jn £ Jn Jn 


where 


<|d| + |/?|- 
£ 


IlFlIc 


0 ‘'O 

T r T 


f ec lf ( aMdr d( + -[ b(s ) [ F(Oe ' f ( aMdr d( 
Jn £ Jn Jn 


ds 


= \d\ + \/3\- 


l|F||, 


■ f 

Jo 




+ 


-Halloo J h(s)[l-e~]d5 


< 


|d| + m l [llFlL ^(l - e=f)] + i IIFIL J b(s ) [l - e^J d. 

= |d| + W\ [||F|L i (l - e=f)] + i IIFIL J b(s ) [l - ef] ds 

<|d| + |jS|||F|L- + -||F|| 00 ||6|| 1 <C 1 

a a 

[ h(s) ds = llftlli- 

J n 


Consequently, we obtain the following inequality: 


u(0) < C 


and 


—ax / —aT \ 

|u(x)| < |u(0)| e — + (1 - e~Ja _1 ||F|L 

< [C 0 _1 |d| + C 0 -fm + WbWJa-'WFlL] [e^ + (l - e=^) otHIFIL] 

< [C 0 _1 |d| + CifHljffl + l|6|| 1 )a -1 l|F|| 00 < C 0 . 

We now prove the estimate the Equation (4): 

\£\ < C, |u(x)| < C 0 . 

If the above inequalities replace the Equation (6), we have the following inequality 

—ax / —aT \ 

|u(x)| < |u(0)| e — + (^1 — e~j a _1 ||F|L (10) 

where |u(0)| < C, a _1 ||F|L ^ C. 

If we take derivative of the Equation (10) we obtain as following: 

COC ~ ax CC ~ ax ( 1 —ax \ 

|u'(x)| < |u(0)|— e~ H— e~ < C 1 H— e ~) 

£ £ \ £ / 

thus, the proving of the Lemma 1 is completed. 

3. The Establishment of Difference Scheme 


Let us consider the following any non-uniform mesh on [0, L\, 

(i) N = {Xq < x t < ■■■ < x w _! < x N , h t = x t — X t _i) 

CO N = (i) N U {x = 0 , X = L}. 

First, we will construct the difference scheme for the Equation (1). First, we integrate the Equation (1) over 

(Xi_i,Xi) 
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*t t r x i _ 

su (x))dx + hf 1 I f{x, u)dx = 0, i = l,N — 1 (11) 

i-l ■'JCj-i 



and we obtain the following inequality 


h t 1 £ | u (x;) 


rXi rX 

^(^i-l)) + ^i _1 I f(x,u)dx = £U fi + hf 1 


f{x,u)dx = 0 


where 


f x ‘ 

R i = ~ h i 1 O; - *i-i)/ (*. 


u)dx. 


So, from the Equation (12) the difference scheme is defined by 


eu fi + fdXi.Ui) + Ri = 0, i = 1,1V 

Second, we define an approximation for the boundary condition of the Equation (1). 

w x 

u 0 = fdu N + / I b(s)u(s)ds + d 

x- N X- 

^ J b(s)u(s)ds = ^ J [b(s)u(s ) — bjiqjds 

• - r ' • i=1 

W « 

biUihi = r + ^ hjUjh; 


where 




i=l 


and remainder term 


JV 

•If 


=1 1 


(x - x^) — [b(t)u(t)]dt. 


u 0 


— fdu N + ^ hibiUi + d + r. 

i =1 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


(17) 


Thus, by neglecting and r in the Equation (13) and the Equation (17), we suggest the following difference 
scheme for approximating the Equations (l)-(2): 

(18) 


£yx,i += o,i = i,N 


y o 


= PVn +'^ j h i b i y i + d 

i=1 


(19) 


4. Bakhvalov Mesh 

We use Bakhvalov mesh, because the difference scheme (18)-(19) must be £ —uniform convergent. 
According to the transition point following as: 

cr = min j—, — a _1 £ln£j 

the piecewise uniform mesh divides each of the interval [0, cr] and [cr, L] into j equidistant subintervals. 
The mesh point x t is defined as follows: 


a < — ise, X; G [0, cr], x ; = —a 1 £ln 

-i. 


cr — — ise , Xj G [0, cr], x t — —a 1 £ln 


N 
’ 2 


t; G [cr,L], Xj = a + (^i — h, i = ^ + 1,..., N, h = 


1 - (1 -e)-]. i = 0,1,...,- 


2(L-o) 
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5. Uniform Error Estimates 

In this section, we obtain the convergence of the method. First, we give the error function 
Zj = y t — U;, i — 0,N, z t is the solution of the discrete problem, 

sz x -,i + fdXi.yt) - fdXi.Ui) =R it i = VAT (20) 

N 

Z 0 = P z n +'Y j h i b i z i -r (21) 

i =1 

where R t and r are defined by Equation (12) and Equation (16). 

Lemma 5.1. Let z, be the solution to Equation (20), 

|b*| = maXjlhjl, p k — '-j-, k — 1-2 and a — max ; |a ; | 


and 


Then the estimate 



llz|loo, Sjv < C(|| R IL )5Jjv + |r|) 


( 22 ) 


holds. 


Proof. We use intermediate value theorem to rewrite the Equation (20) as follows: 

-/Ui.^i) = /u(yi-Ui) = f u Zi 

Here we assumed that 

fu ~ a i- 

Thus, we can write 

ez x ; + a^Zi = R u i = 1,N — 1 

(23) 


N 

z 0 = /3z n +^h i b i z i -r. 

i=1 

(24) 

Rearranging Equation (23) gives 

and 

e h i R i 

z, ~ - T z i- 1 1-r~ 

£ + a^i £ + a^i 



z i =z 0 [] _ |<?fc]+^%[] — [ Rj] 
k= 1 k= 1 j=k+ 1 


where 

£ 

q k = -— 

£ + CL^ rife 



and 


h k R k 

(p k =-—. 

£ + CL k tl k 

According to the maximum principle for Equations (23)-(24), we have the following inequality 

\\ z \L,a> N < |z 0 1 + a -1 ||f?lloo,s w ■ 
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Now, we evaluate the Equation (24). After a simple calculation we obtain 

N 

\z 0 \ < p\z N \ + + |r| (26) 

i=1 


where 


\b t \ = max|bj| 


and 


N 

Z N = Z O 0 N +^<p k 0 N - k . 
k=1 


The Equation (27) write in the Equation (26), we have 


Zol ^ P 

N 

Z O 0 N + S ^_ i <Pk®N-k 

N 

+ b t '^ j h i 

i 

z o 0 i +^ j cp k 0 i -k 


k= 1 

i= 1 

k = 1 


+ r 


After doing some arragement, 


\z 0 \ < 


P\Yf k=1 (Pk®N-k \ + b,Yf =1 K\Y! k= i (p k 0j- k \ + 
(l-/?|0 w |-fc t 2^10,1) 


it follows that 


kol <■ 


St- tst) (n^)“ I + »■ S'- h ‘ a-. feS:) fck) ' 


1+apkJ \l+ap k > 


+ r 


< CdlKlL,^ + |r|) 


where 


p k = —, k — 1-2 and a = maXjlaJ. 

If the Equation(29) write in the Equation (25), we have 

llzlloo,^ < C(l|f?lloc,6j w + M) + « -1 |l^lloc,6j w 

< COlKlUs* + |r|). 

Lemma 5. 2. Under the assumptions of Lemma 1 and j- > a > 0, f(x,u ) 6 C 2 [0,L], the solution of (l)-(2) 
satisfies the following estimates: 

II R\Lfi M < CN- 1 (30) 

|r| < CN- 1 (31) 

Proof. First, we consider error function R h 


(27) 


(28) 


(29) 


r*i 

|/?J < hj 1 I (x; — X;_i)|/ (x,u)|dx 

r x i 

= h i 1 I (x f - Xi_i) 

J Xi -1 


df dx df{x, u(x)) du 
du dx du dx 


dx 


r x i r x i 

< Ch.- 1 I (x; — X; _f)dx + Chj" 1 I (x; — Xi-i)^' 

‘'r i dy. 1 


(x)|dx 
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1 Zf£ 

11 —e £ 
£ 


r ‘ 

< Chi + Ch t 1 (x t — Xj-i) 

r _ r -ax ] 

| h i +h i 1 £~ 1 I (Xj — x i _ 1 )e £ dx), 


dx 


< C 


1 <i<N. 


(32) 


In the first case, x ; 6 [0, cr], 


a) cr < cr = —a 1 sln£: 


hj — X; — x^! = a x £ 


In 


1 - (1 - £ ) 


2 i 
N 


-In 


1 - (1 - £ ) 


2(t~ 1) 
N 


where, according to t, we use intermediate value theorem as follows: 

2(1 - £)iV _1 

hj = a -1 £--— 7 ;- . — 1 < 2a _1 (l - £)iV -1 < CN _1 

1 1 - ii2(l - e)^- 1 

Thus, we obtain 

{ C x i —ax 'j ( r x i —ax 'j 

hi + hj“ 1 £ _1 I (x; — x i _ 1 )e~dx) < C jhj + h i “ 1 £ _1 h i I e~dx) 

< C [hi + e" 1 J 1 e^dxj = C (hi + cT 1 [e - ^ - e - ^]] 


(33) 


e a 1 £ln[t-(l-£)| L ] _ g a 1 £ln[l-(l-£)^ r 1) 


= C |dj + a x e £ 

< CiCN- 1 + 2a _1 (l - s^N- 1 } < CN- 1 . 


]]} = C{hi 


+ 2cr 1 (l-£)W" 1 } 


b) cr — - ,ct — - < —a ^^slns : 
J 2 2 


h L — X; — Xj_ x = a 1 s 


In 


i _( i _ e2 :)|]_ in [ i _( i _ e2 :) 


=£a 2 (i - 1) 


IV 


where, according to i, we use intermediate value theorem as follows: 


2(1 - e^)N~ 1 

v _ 

l-t 1 2(l-e^)iV- 1 


^ = a^e --- : —^-< 2a- 1 N- 1 < CN- 1 . (34) 


Hence 


(r ~ ux j ~ ux j- i ~n l ^ 

hi + a -1 |e £ — e £ Jj = C jd; + a _ 1 e~ 
< CAT 1 

In the second case, x ; G [cr, 1], we obtain 


12 Mn|l-^l-e2£j-j « £ln|l-(l-e 2 £ j N 


\R t | < C j hi + E- 1 j 1 eTTdxj < C + e" 1 ^ j < CfAT 1 + e^/T 1 } < C{N- 1 + e” 1 AT 1 } < CN- 1 . 


Consquently, we have 

II fllL,^ < CN- 1 . 

Now, we evaluate error functionr as follows: 
x — Xi_ x < h it b* — max|b(x)| < C , b (x) < C ve u (x) < C 0 
under the above assumptions, 


(35) 
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r < 


N 

If 

i=1 x i~ 1 


(x — x^) [b'(x)u(x) + b(x)u (x)]dx 


N 

i—1 x i~ 1 
N 

<c£>,r 

U V, 


1 

1 + -e 
£ 


e £ 


dx 


dx 


< C 


5 >f 

1 1=1 x i~ 1 


1 

1 + -e 
£ 


dx + 


r 


A 


1 H—e f 

dx 

j x , 

a i-i 

£ 

y 


In the first case, x,- 6 [0, cr]: 
a) cr < -, cr = —arsine. 


[ r ff r 1 <«] 

r 1 r l “1 


d; 1 + - e £ 

dx + hi I 1 + — e £ 

dx 

[ Jo L £ 

L £ 



< Ch , 


(36) 


W-l 


- 


a 1 £ln 


1 - (1 - £) 


2 (i - 1) 


IV 


and 


hj — X; — Xj_! = a 1 £ 


In 


1 - (1 - £) 


2i 

iv 


-In 


1 - (1 - £) 


2(t~l) 

N 


where, according to i, we use intermediate value theorem as follows: 

1 2(1 — £)iV 1 


h[ = a 1 £- - ) — < 2a *(1 - e)1V 1 < CIV 1 

1 1 - ^ 2(1 - s^N- 1 


Thus, we obtain 


|r| < C[h t + < CN' 1 . 


1 1 

b)cr = — ,cr = — < —a l £ln£, 

J 2 2 

After a simple calculation we obtain. 


. 2(l-e^)iV- 1 , 1 

^ = a _1 £----< 2cc~ 1 N~ 1 < CN- 1 . 

1 - ii2(1 - e^N- 1 


Hence 


|r| < C[hi +h t ] < CN 


-l 


In the second case, x t 6 [cr, 1]: 

. 1 _1 , 1 
a) cr<-,—a 1 £in£ <-. 
' 2 2 

we use the following equation. 


, 2(1-cr) 2 N 

h = hi = h 2 — -—-< —, i = — + 1, N 

1 N N 2 


we obtain. 


(37) 


(38) 


|r| < C[hi + h^ < Cih 1 + h 2 ] < C^” 1 + 2 AT 1 ] < CiV" 1 
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\r\ < CN 


-l 


(39) 


Teorem 5.1. Under the assumptions of Lemma 5.1. and Lemma 5.2. Let u(x), be the solution of (l)-(2) and y h 
be the solution of (18)-(19). Then, the following uniform error estimate satisfies 


lly- wIIoo.sa, <cn 1 . 


(40) 


6. Algorithm and Numerical Results 

We give some numerical results to solve the nonlinear problem (18)-(19) on Bakhvalov mesh using 
quasilinearization technique. 

6.1. Algorithm 

We present the following quasilinearization technique for the difference scheme (18)-( 19): 




Cn—1) ( 




i = 1 


Vi 


(n) _ 




_ (£. , d/Oj.yj) l 

\h i dy ) 


i = 1 ,..,N 


yo n) = PVn 1] + > i; + d, n = 1,2,... 


N 


(n—1) 


i=l 


6.2. Numerical Results 

We give some examples to see the effectiveness of the presented method 

Example 1 

eu + 2u — e~ u + x 2 = 0, 0 < x < 1 
1 f 1 


1 If 1 

u(0) = -u(l) + — I e S u(s)ds + 1 

2 4 Jr, 


(41) 

(42) 


this problem has not the exact solution. Thus, we use co —mesh to calculate the errors and convergence rates, 
respectively. 


e N = max ef , e" = max|y" - y 2 2 "|, P E N = log 2 


where, y- v and yf/ v are the numerical solutions for IVand 2 N. 
According to (6.1), we give the algorithm as follows: 


^ ( n ) (^ I o _i_ —y^ n M n 0 

-{h, +2+e ■ ) y ‘ =2y 


, 2 (n-1) . -y. Cn_1) '\ -y^-V 

+ x~ — y ; I 2 + e y <- I — e 


y /" 3 = - 


£ (n) (n- 1) . 2 (n-1) 

y,-i - (2y> +xf-yl [2 + e y i J - e ^ ) 


- + 2 + e“ y ; 

hi 


(n-1) 


, i = 1,.., TV 
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N 

y 0 (n) = ^y« n_1) + X +!’ n = !' 2 ' - 

i =1 

the initial guess in the iteration procedure is y® = 0.5, 

The stopping criterion is taken as |— yf 1 ' ) | < 10 -5 . 

The numerical results obtained from the problem of the difference scheme by comparison, the error and 
uniform rates of convergence were found and these are shown in Table 1. Consequently, numerical results show 
that the proposed scheme is working very well. 


Table 1. The computed maximum pointwise errors e N and e 1N , the numerical rate of convergence p N on the 
Bakhvalow mesh co N for different values of Wand £. 


s/N 

8 

16 

32 

64 

2- 1 

0.0364705764 

p=0.76 

0.0213939341 

p=0.88 

0.0115893393 

p=1.05 

0.0051560120 

2~ 2 

0.0457353908 

p=0.79 

0.0263063297 

p=0.89 

0.0141504831 

p=0.94 

0.0073451227 

2~ 3 

0.0547795035 

p=0.79 

0.0316355902 

p=0.87 

0.0172632555 

p=0.93 

0.0090159047 

2_ 4 

0.0592489265 

p=0.74 

0.0354285745 

p=0.87 

0.0192899598 

p=0.93 

0.0100776157 

2 -5 

0.0617440905 

p=0.71 

0.0376291853 

p=0.87 

0.0204860025 

p=0.93 

0.0107447933 

2 -6 

0.0632177961 

p=0.70 

0.0389122636 

p=0.87 

0.0212015157 

p=0.92 

0.0111341923 

2~ 7 

0.0640959677 

p=0.69 

0.0396529864 

p=0.87 

0.0216490391 

p=0.93 

0.0113563269 

2 —8 

0.0646130925 

p=0.68 

0.0400746536 

P=0.87 

0.0219009225 

p= 1.00 

0.010926500 



x 


Figure 1. Exact solution distribution for N = 16 and s = 2 4 ,2 6 ,2 8 
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Figure 2. Error distribution lor JV = 16 and e= 2 2 ,2 4 ,2 6 ,2 8 

From the graps it is show that the error is maximum near the boundary layer and it is almost zero in 
outer region in the Figure 2. When s — values are small, graphic curves leaned more towards the coordinate axes 
in Figure 1. 

Consequently, the purpose of this study was to give uniform finite difference method for numerical 
solution of nonlinear singularly perturbed problem with nonlocal boundary conditions. The numerical method 
was constructed on Bakhvalov mesh. The method was pointed out to be convergent, uniformly in the e- 
parameter, of first order in the discrete maximum norm. The numerical example illustrated in practice the result 
of convergence proved theoretically. 
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